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ABSTRACT 

The most important cooling lines of the neutral interstellar medium 
(ISM) lie in the far-infrared (FIR). We present measurements by the Infrared 
Space Observatory Long Wavelength Spectrometer of seven lines from neutral 
and ionized ISM of 60 normal, star-forming galaxies. The galaxy sample 
spans a range in properties such as morphology, FIR colors (indicating dust 
temperature), and FIR/Blue ratios (indicating star-formation activity and 
optical depth). 

In two-thirds of the galaxies in this sample, the [C II] line flux is proportional 
to FIR dust continuum. The other one-third show a smooth decline in L[cii]/Lfir 
with increasing Fy(60 fj,m)/F^{100 /im) and Lfir/Lb, spanning a range of a 
factor of more than 50. Two galaxies, at the warm and active extreme of the 
range have L[cii]/Lfir < 2 x 10~^(3(J upper limit). This is due to increased 
positive grain charge in the warmer and more active galaxies, which leads to less 
efficient heating by photoelectrons from dust grains. 

The ratio of the two principal photodissociation region (PDR) cooling lines 
L[oi]/L[cii] shows a tight correlation with Fy(60 /im)/Fj,(100 /xm), indicating 
that both gas and dust temperatures increase together. We derive a theoretical 
scaling between [N II] (122 /im) and [C II] from ionized gas and use it to 
separate [C II] emission from neutral PDRs and ionized gas. Comparison of 
PDR models of Kaufman et al. (1999) with observed ratios of (a) L[oi]/L[cii] 
and (L[cii] + L[oi])/Lfir and (b) L[oi]/Lfir and F^(60 /im)/Fj,(100 fim) yields 
far-UV flux Go and gas density n. The Gq and n values estimated from the two 
methods agree to better than a factor of 2 and 1.5 respectively in more than 
half the sources. 

The derived Gq and n correlate with each other, and Gq increases with 
n as Gq (X n'^, where a ~ 1.4 . We interpret this correlation as arising 
from Stromgren sphere scalings if much of the line and continuum luminosity 
arises near star-forming regions. The high values of PDR surface temperature 
(270 — 900 K) and pressure (6 x 10^ — 1.5 x 10'^ K cm^^) derived also support the 
view that a signiflcant part of grain and gas heating in the galaxies occurs very 



- 2 - 



close to star-forming regions. The differences in Gq and n from galaxy to galaxy 
may be due to differences in the physical properties of the star-forming clouds. 
Galaxies with higher Gq and n have larger and/or denser star- forming clouds. 

Subject headings: radiation mechanisms: thermal, ISM: atoms, ISM: general, 
ISM: H II regions, galaxies: ISM, infrared: ISM: lines and bands 

1. Introduction 

The atomic and ionic fine-structure lines [C II] (158 fim) and [O I] (63/im) are the 
dominant cooling lines for neutral interstellar gas, and ionic fine-structure lines like [O III] 
(88 /xm), (52 /xm) and [N II] (122 /xm) are strong coolants in H II regions. These lines can be 
used as diagnostics to infer physical conditions in the gas, such as temperatures, densities 
and radiation fields, by comparing with models of photodissociation regions (PDRs e.g. 
Tielens & HoUenbach 1985, Sternberg & Dalgarno 1989, Wolfire, Tielens & HoUenbach 
1990, Kaufman et al. 1999) and H II regions (Rubin et al. 1991). 

Among these, [C II] is the most ubiquitous and best studied line. It was predicted to 
be the dominant coolant for diffuse neutral media by Dalgarno & McCray (1972), but the 
first detection of this fine was towards the dense star-forming regions M42 and NGC 2024 
(Russell et al. 1980). The early observations of external galaxies were of nearby or IR 
bright galaxies and showed that [C II] emission was associated with dense gas irradiated by 
ultra-violet (UV) light from young star-forming regions, often near galactic nuclei (Crawford 
et al. 1985, Stacey et al. 1991). Later, observations of quiescent galaxies like NGC 6946 
showed that averaged over the whole disk, including the Hl-rich outer disk, a significant 
fraction of the total [C II] in a galaxy may arise from diffuse ionized or diffuse atomic gas 
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(Madden et al. 1993, but see Contursi et al. 2001). Observations of many FIR lines in the 
Milky Way by FIRAS , including [C II] (158 /im) and [N II] (205 fim), indicated that a fair, 
and possibly dominant, fraction of [C II] emission is from diffuse {rif, = 1 — 5 cm~'^) ionized 
gas (Petuchowski & Bennett 1993, Heiles 1994). Recent ISO observations have yielded 
more surprises. In a previous paper we reported a deficiency of [C II] compared with FIR 
continuum in three normal star-forming galaxies (Malhotra et al. 1997, Paper I). A similar 
deficiency was reported for ultraluminous galaxies by Luhman et al. (1998). 

The range of interpretations for the origin and behavior of the [C II] line in particular, 
and FIR cooling lines in general, emphasizes the need for studying a suite of FIR lines in a 
large sample of galaxies. The Long Wavelength Spectrometer (LWS, Clegg et al. 1996) on 
the Infrared Space Observatory (ISO, Kessler et al. 1996) has made it possible to observe 
a large number of atomic and ionic fine-structure lines in the FIR with unprecedented 
sensitivity, so that a large sample of galaxies could be observed. 

In this paper we report on and interpret observations of fine structure atomic and ionic 
lines observed in 60 normal galaxies. The sample consists of galaxies whose luminosity is 
dominated by star- formation and excludes Active Galactic Nuclei (AGN). We observed 
60 distant galaxies, for which all the FIR emission is within one LWS beam, as well as 6 
nearby, resolved galaxies. Here we report on the distant sample. 

We observed from one to seven atomic and ionic fine structure lines for each 
galaxy depending on its FIR brightness. The fines include [C Il](158//m), [O l](145//m), 
[N Il](122//m), [O III] {88fxm), [O l](63//m), N III (57//m) and [O III] (52//m). The most 
common lines observed are [C Il](158/im), [O I](63/xm), [O III] (88;Lim) and [N Il](122/xm). 
There is no convincing detection of N III (57;um), so we do not discuss that line here. This 
paper offers a look at the statistical behavior of [C II], [N II], [O III] and [O I] lines in a 
diverse sample of galaxies, aimed at better understanding the physical conditions in the 
Interstellar medium (ISM) of normal star-forming galaxies. 

The paper is arranged as follows: in Section 2, we discuss the far-infrared lines 
observed, their properties, and their diagnostic value in deriving the physical conditions 
in the ISM. In Section 3 we discuss the observations and data reduction procedures. In 
Section 4 we describe the statistical trends observed; in Section 5 we interpret these trends, 
and in Section 6 we discuss the physical conditions in the PDRs in these galaxies. Section 7 
contains conclusions and a summary of the main results. Appendix A contains a description 
of the sample selection and how the sample spans the parameter space in galaxy properties. 
Appendices B and C contain tables of line fiuxes and derived physical quantities: far-UV 
fiux Go, gas density n, temperature T and pressure P. 
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2. Fcir Infrcired Fine structure lines 
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Properties of FIR lines observed 
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NIII 


57.317 


29.60 


47.45 




3 X 10=^ [e] 1 





" Potential required to create the ion 
Potential required to ionize the species 
critical density 

1 Genzel 1991 

2 Hollenbach & McKcc 1989 

3 Launay & Rocff 1977 
^ Jaquet et al. 1992 

Table 1 summarizes the properties of the lines observed for this sample and discussed 
in this paper (cf. Genzel 1991 for a more exhaustive description of atomic, ionic and 
molecular lines in the infrared) . Here critical density for a transition is the density at which 
collisions balance spontaneous radiative transitions. T — AE/k indicates the temperature 
corresponding to the energy difference in the upper and lower levels. 

2.1. [C II] 

The C+ fine structure transition at 157.714 //m is the most important coolant of the 
warm neutral interstellar medium. Carbon is the fourth most abundant element and has 
a lower ionization potential (11.26 eV) than hydrogen, so that carbon will be in the form 
of C+ in the neutral surface layers of far-UV illuminated neutral gas clouds. The depth of 
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these C+ zones is generally determined by dust extinction and often extends to Ay < 4. 
The 158 /xm [C II] line is also relatively easy to excite (AE/k ~ 91 K), so that C+ can cool 
warm {30 K < T < 10"^ K) neutral gas where the two most abundant atoms H and He 
cannot (cf. Tielens & Hollenbach (1985), hereafter TH85; Wolfire, Tielens & Hollenbach 
1990, hereafter WTH90). In the zones of PDRs, the electrons are supplied by the 

ions so that the abundance of electrons relative to hydrogen is x{e) ~ 10~^. In such 
conditions the [C II] (158 fim) line is excited by coUisions with hydrogen atoms 

[C II] from neutral regions can be used in combination with the [O I] lines and the 
FIR continuum, to derive the gas density n and the incident FUV (6 eV < hu < 13.6 eV) 
radiation flux Gq respectively (e.g. Wolfire, Tielens & Hollenbach 1990; Kaufman et al. 
1999). Go conventionally is the FUV flux normalized to the average local interstellar flux of 
1.6 X 10"^ ergs cm~^ s~^ (Habing 1968). Since [C II] comes from both ionized and neutral 
regions, decomposing the two components may be necessary before the comparison with 
[O I] can be made. 

2.2. [N II] 

Since nitrogen has an ionization potential (14.5 eV) higher than that of hydrogen, 
[N II] (122 iim) arises only in ionized gas. This flne structure fine is excited by coUisions 
with electrons. [N II] (122 iim) has a critical electron density of 3.1 x lO^cm"^, and [C II] 
(158/im) has a critical electron density of 50 cm~^ (Table 1). We can estimate the [C II] 
arising from low density (rif. < 50cm~^) ionized regions by comparison with the [N II] (122 
/im) line, as described in section 5.4 below. 

2.3. [O I] 

Oxygen has an ionization potential of 13.62 eV, quite close to that of hydrogen, so 
atomic O is found in neutral regions only. In PDRs 01 can exist in atomic form far deeper 
into clouds than C+. All oxygen not incorporated into CO can stay atomic to depths as 
large as Ay = 10 if relatively high FUV fluxes impinge on a cloud (TH85). 01 has two 
flne structure transitions, at 63/im and 145;um. The excitation energy of [O I] (63 fim) 
corresponds to 228 K and the critical density is ~ 5 x lO^cm"^ (at T ^ 300 K), so 01 
lines arise in warm and dense neutral regions. As discussed above, [O I] (63 /xm), [C II] (158 
//m) and the FIR continuum constrain Gq and n. The 145//m line lies AE/k — 325 K 
above the ground state so the ratio of [O l](145/xm)/[0 I] (63 /xm) measures temperature 
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(or Go) in the temperature range of ~ 300K. In many cases, however, this ratio may be 
an indicator of optical depth in the 63;um hne, which is often moderately optically thick. 
Early observations and modeling of the [O I] lines from Orion indicated that the 63 //m line 
was optically thick in emission from the Orion PDR (Stacey et al. 1983, 1993), while more 
recent observations of Sgr B2, DR21 and NGC6334 indicate self absorption in the [O I] 
line by intervening cool O along the line of sight (Keene et al. 1999; Poglitsch et al 1996; 
Kraemer et al. 1996). 

2.4. [O III] 

Since O II has an ionization potential of about 35 eV (Table 1), O III lines come 
from ionized regions of galaxies where the UV radiation field is dominated by fairly early 
type (~06) stars. The ratio of the two lines of [O III] at 88 /im and 52 /im is useful 
for determining electron density in H II regions. The rie derived is insensitive to the 
temperature of electrons Tg (Rubin et al. 1994) and independent of abundance variations. 
Since the line ratios give a measure of local density, a large beam measurement by this 
technique gives < rig >, as opposed to < >^^^ provided by emission measures and 
observed extent. 



2.5. NIII 

The ratio N IIl(57//m)/[N II] (122 //m) provides a measure of the effective temperature 
of the ionizing star(s), Tg// (Rubin et al. 1994). N IIl(57//m) is only marginally detected in 
some sources and we are not sufficiently confident of the measured fiuxes to include N III in 
this paper. 

3. Observations and data analysis 

All the line observations for this project were obtained using the low-resolution grating 
mode of LWS (Clegg et al. 1996). With the LWS spectral resolution of 0.29 and 0.6//m 
(for wavelength ranges 43-90.5 /im and 90.5-197 fim, respectively) we do not resolve the 
lines. Sixty distant galaxies were selected to have FWHM < 30" in FIR emission using 
deconvolved IRAS maps . Therefore, with the LWS 70" beam, we report the total line 
fiuxes for these galaxies, and, using measured distances, the total luminosities from these 
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galaxies 0. The observations were made in LWS02 mode where the gratings are scanned 
across known hnes. A few galaxies were observed in the LWOl mode where a complete scan 
was made to cover the wavelength range 45-193 fim. 

The hne observations were planned to achieve a sensitivity of a fixed fraction of the 
FIR continuum flux, depending on the expected line strength and the feasibility of the 
observations. Between one and seven of the lines [C II] (158 /xm), [O I] (145 yum), [N II] (122 
/xm), [O III] (88 /xm), [O I] (63 /im), N III (57 fim) and [O III] (52 fim) were observed for each 
of the 60 galaxies in the sample, depending on what was practical given the FIR brightness 
of the galaxy. The most commonly observed lines were [C Il](158/im), [O l](63/im), [O III] 
(88/im) and [N II](122/xm), in that order. Table 2 gives the planned sensitivities for the 
different lines along with the 1-cr noise actually achieved for these lines averaged over the 
whole sample of galaxies. 



Table 2: Planned and achieved line sensitivities 



Line 


wavelength 
/im 


l-cr flux (planned) Average 1-cr flux (achieved) 
(as a fraction of the continuum flux, Ffir) 


[C II] 


158 


2.2 X 10"^ 


1.4 X 10-^ 


[0 I] 


145 


2.0 X 10~^ 


3.0 X 10"^ 


[N II] 


122 


8 X 10-^ 


7.0 X 10-^ 


[0 III] 


88 


1 X 10-^ 


1.3 X 10-^ 


[0 I] 


63 


1.0 X 10"^ 


1.3 X 10-^ 


NIII 


57 


1.3 X 10"^ 


2.1 X 10^^ 


[0 III] 


52 


2.4 X 10"^ 


2.4 X 10^^ 



The [C II] line observations were planned to achieve (la) sensitivities of 2.2 x 10 ^ x Ffir, 
where FpiR is the far-infrared flux of the galaxy between 42/im and 122/im and is computed 
according to the relation FpiR = 1.26 x 10-^^[2.58 x F^(60/im) + F^{lOOfxm)]Wm-'^ 
(Helou et al. 1988), where F^{60fim) and -Fjy(100/im) are flux densities in Jansky in the 
IRAS 60 and 100 /im bands. For comparison, previous observations of the Milky Way, 
starburst galaxies and galactic nuclei show that the line to continuum luminosity ratio 
L[cii]/Lfir = 1 - 10 X 10^3 (Stacey et al. 1985, Wright et al. 1991, Crawford et al. 1985, 
and Stacey et al. 1991). In Galactic PDRs associated with H II regions L[cii]/Lfir varies 
from ~ 3 X 10~^ in NGC 2023 (a reflection nebula) to ~ 8 x 10"^ in W51, decreasing with 
H II region or PDR density and FUV flux (cf. Crawford et al. 1985, HoUenbach, Takahashi 
& Tielens 1991). 



In practice we will always use ratios, so the distance errors cancel 
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The data were calibrated with the ISO pipehne OLP7.0 and reduced by two methods. 
In the first method the hne profiles were derived from several scans by running a median 
boxcar filter through them. An automatic rejection algorithm was used which discarded 
measurements more than 5cr away from the mean flux level and further discarded four 
measurements subsequent to the discarded points. This sigma-clipping gets rid of cosmic 
rays and a time-dependent rejection is done to remove possible memory effects on the 
detectors. The median of the observed fluxes is used instead of the mean to reduce 
the influence of outlying points arising from cosmic ray hits. The flux in the lines was 
determined by directly integrating under the line, after fitting a linear baseline to the 
underlying continuum measurements. The upper limits on non-detections were derived by 
calculating the flux from a hypothetical gaussian line with an amplitude of 3a and the 
effective instrumental profile, since the lines are unresolved for all sources. In the second 
method, the individual scans were inspected and bad points were rejected after manual 
inspection. A linear baseline and a gaussian profile was then fit to the data. The line fiuxes 
are derived from the fitted gaussian. The fiuxes from the two methods agreed for well 
detected (> 5a) sources. The hne fluxes thus derived are tabulated in Appendix B. 

Presently, the calibration uncertainties in the fiux measurements are due to two sources: 
(1) ill-determined dark currents, which are additive in nature and affect the measured 
line fluxes minimally (< 5% in most cases) because the dark currents are subtracted with 
the continuum levels while determining the line flux; and (2) the absolute calibration 
and spectral response is tied to observations and a model atmosphere of Uranus. The 
relative flux calibration is determined by a measurement of a few bright lamps in between 
observations. We estimate that the fiux uncertainty in our sources is better than 20%, 
including all random errors and systematic calibration terms. The instrument team reports 
calibration better that 10-15% for compact, bright sources observed on-axis (Swinyard et 
al. 1999). 

4. Observations and Phenomenology 

4.1. The [C II] (158 /im) observations 

While [C II] and [O I] are the main cooling lines of the neutral atomic ISM, the gas 
heating is dominated by photoelectrons from dust grains (Watson 1972; for a recent review 
see HoUenbach & Tielens 1999). In a fairly indirect and inefficient mechanism, incident 
FUV photons with energies high enough to eject electrons from dust grains (h. u > 6 eV) 
heat the gas via these photoelectrons, with a typical efficiency of 0.1 — 1%. Efficiency is 
defined as the energy input to the gas divided by the total energy of the FUV photons 
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absorbed by dust grains. This efficiency is determined by the microphysics of the grains, 
in particular the work function, photoelectric yield and the charge of the grains (which is 
determined by the ratio of FUV fluxes to gas density Go/n). Before ISO observations little 
variation was seen in the total L[cii]/Lfir of galaxies (Stacey et al. 1991, Crawford et al. 
1985), or in the Milky Way at 7° scales (Wright et al. 1991). This suggested that averaged 
over galaxy scales, the heating efficiency did not vary much, and guided our line sensitivity 
choices, which were scaled to the FIR continuum. 

The two circumstances where we do expect the [C II] /FIR ratio to change are: (1) 
for a high ratio of FUV flux to gas density, Go/n, the grains get positively charged raising 
the potential barrier for photoelectric ejection, thereby dropping the heating efficiency; and 
(2) when the hardness of the radiation changes, changing the ratio of FUV light which is 
effective in photoelectric heating the gas to less energetic light which can only heat the 
dust. In subsequent sections we will see examples of each of these. 

In a previous paper (Paper I) we presented [G II] measurements for half the current 
sample, and saw large variations in L[cii]/Lfir including three non-detections of the [G II] 
fine in sources with relatively large FIR continuum fluxes. With a full sample of 60 galaxies 
and with better calibrations we verify the main observational results of Paper I. From our 
comprehensive data (Table 6 and Figure 1) we observe: 

(1) About two-thirds of the observed galaxies (41 out of 60) show a ratio 
L[cii]/Lfir > 0.2%. This is consistent with models of photoelectric heating in PDRs 
illuminated by moderate FUV fluxes (TH85) and with previous observations (Stacey et al 
1991). 

(2) There is a trend of decreasing L[cii]/Lfir with warmer FIR colors, 
Fi,(60 //m)/Fi,(100 /im) and increasing star-formation activity, indicated by higher 
Lfir/Lb ratios, where Lb is the luminosity in the B band. Wc test the significance of 
these correlations by performing the generalized Kendall's test (Isobe, Feigelson & Nelson 
1986, Brown, Hollander & Korwar 1974) which accounts for censoring of data (i.e. upper 
limits). Since we don't know the exact functional form of the dependence of L[cii]/Lfir 
on F,y(60 /xm)/F,y(100 /^m) and Lfir/Lb, a rank test is appropriate. That L[cii]/Lfir and 
F^(60 /xm)/F^(100 /xm) are uncorrelated is excluded at the 99.999% (4.4 a) level. The 
anti-correlation between L[cii]/Lfir and Lfir/Lb is somewhat weaker and the hypothesis 
that L[cii]/Lfir and Lfir/Lb are uncorrelated is rejected at the 99.5% level (2.8(t). 

(3) Three galaxies (NGG 4418, IG 0860 and GGGG 1510.8+0725) near the extreme 
end of this trend in the Fj^(60 /xm)/Fj^(100 /xm) ratio and in the ratio of Lfir/Lb showed 
no detectable fine emission in the [G II] and [O I] fines, as reported in Paper I. Here we 
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report the observation of CGCG 1510.8+0725 with greater sensitivity and the detection 
of the [C II] hne, consistent with the 3a upper hmit reported earher. No other prominent 
hnes, e.g. [O I] (63 //m), are seen in these galaxies. 




0.4 0.6 0.8 1 1.2 1.4 -1 

F,(60 ^m)/F,(100 ^m) Log(Lp,j,/LB,^J 



Fig. 1. — (a) The ratio of [C II] to far-infrared continuum, L[cii]/Lfir, is plotted against the 
ratio of fiux in the IRAS 60/xm and 100/im bands, F,^(60 /xm)/F,^(100 fim). Filled circles 
are normal galaxies from the ISO-KP sample and triangles are luminous and ultraluminous 
galaxies from the sample of Luhman ct al. (1998). The line fluxes arc uncertain by about 
20%. There is a trend for galaxies with higher F,^(60 /.im)/Fj^(100 /xm) (indicating warmer 
dust) to have lower Ljcuj/LpiR, for normal as well as ULIRGs. Two normal galaxies in a 
sample of 60 have no detected [C II] , and they are identified with labels and shown as upper 
limit symbols in the figure; other upper limits come from Luhman et al. (1998). Rank 
correlation tests show that L[cii]/Lfir and Fi,{QO /im)/F,y(100 /xm) are correlated at the 
4.4(7 level, (b) L[cn]/LFiR shows a similar but weaker trend {2.8a signiflcant) with the ratio 
Lfir/Lb, which is an indicator of star formation activity and optical depth in dust. Galaxies 
with higher Lfir/Lb (more active star formation) have lower L[cii]/Lfir. The decrease in 
L[cn]/LFiR with Lfir/Lr is not continuous but sets in for galaxies with Lfir/Lr > 0.8. 

Figure 1 also displays data from the sample of luminous and ultraluminous infrared 
galaxies (ULIRGs) of Luhman et al. (1998). They follow the same trends as normal 
galaxies. 

We also explore the dependence of L[cii]/Lfir on the total FIR luminosity of the 
galaxies (Figure 2). To extend the range of luminosity explored, we plot the normal galaxies 
as well as the ULIRGs from the sample of Luhman et al. (1998). There seems to be a 
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Fig. 2. — The ratio L[cii]/Lfir plotted against the far-infrared luminosity of galaxies, LpiR. 
Filled circles represent normal galaxies from the ISO Key Project sample and triangles denote 
luminous and ultraluminous galaxies from the sample of Luhman et al. (1998). There is 
a trend for galaxies with higher luminosity to have lower L[cii]/Lfir. Two normal galaxies 
in a sample of 60 and three ULIRGs have no detected [C II]; they are identified and shown 
as upper limit symbols in the figure. The correlation between luminosity and Ljcuj/LpiR 
is weakened by having the [C II] deficient normal galaxies in the middle of the luminosity 
range. 

dependence on the FIR luminosity, but many of the [C II] deficient galaxies lie in the middle 
of the luminosity range and the spread in L[cii]/Lfir is large at the high luminosity end. 
With only the sample of normal galaxies the correlation between Lpiq/LpiR and luminosity 
is significant at the 2.1a level; adding the luminous and ultraluminous galaxies from the 
sample of Luhman et al. makes the correlation stronger (3.8 a). We suspect that the 
dependence on luminosity is a secondary correlation and is due to the correlation between 
FIR colors and luminosity (Figure 13b in Appendix A). 

4.2. The [N II] (122 /im) observations 

[N II] at 122 fim arises from ionized gas, both in dense H II regions and diffuse 
ionized gas. Figure 3a shows the decrease of L[nii]/Lfir with increasingly warm FIR colors. 
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Unfortunately the detections are fewer than for the [C II] (158 //m) and [O l](63//m) hnes. 

We detected this hne in about 36% (19 of 52) of the galaxies observed, and the entire range 
of L[nii]/Lfir, including upper limits, is only about a factor of 12. Figure 3a suggests that 
[Nil] follows similar trends as [C II]; i.e. L[nii]/Lfir decreases with F,^(60 /im)/Fi^(100 /xm). 
The anticorrelation between L[nii]/Lfir and F,y(60 /Um)/F,y(100 fim) is 2.5a significant, i.e. 
there is a 1% probability that these two quantities are uncorr elated. 

The ratio L[cii]/L[nii] is fairly constant across a range of FIR colors but with a large 
scatter: L[cii]/L[nii] = 9tg (Figure 3b). A proportionality between the two lines was also 
found in spatially resolved observations of the Milky Way (Bennett et al. 1994). We find 
an average ratio of L[cii]/L[nii]=9 in sources with [Nil] detections, which agrees very well 
with the average empirical ratio in our galaxy (Wright et al. 1991, Bennett et al. 1994). 
However the many upper limits on [Nil] indicate that many sources have larger L[cii]/L[nii]. 




0.4 0.6 0.8 1 1.2 1.4 0.4 0.6 0.8 1 

F,(60 ^m)/F^(100 ^m) F,(60 Mm)/F„(100 ^m) 



Fig. 3. — (a) The [Nil] (122/im) line shows a fairly similar behavior to the [C II] hne, i.e. 
L[nii]/Lfir decreases at high Fi,(60 /xm)/Fi,(100 /xm). This anticorrelation is 2.5(7 significant. 
If true, this suggests either that (much of) [C II] and [Nil] arise from H II regions, or that the 
PDRs which produce [C II] are associated with or surround the ionized gas which produces 
[Nil]. The upper limits shown in this figure are 3a upper limits, (b) This is further illustrated 
by the lack of any observed trends in L[nii]/L[cii] vs F^{60 /xm)/F,^(100 /im). The observed 
L[Nii]/L[cii] = Slg- 
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4.3. The [O I] (63 //m and 145 //m) Observations 

The [O I] line at 63 fim is the second most commonly detected line in the present 
sample. It is detected in 46 of the 53 galaxies where it was observed. It is not detected 
in the [C II] deficient galaxies. Apart from these galaxies, L[oi]/Lfir is unchanging with 
Fi.(60 /im)/F;^(100 /xm). The mean value of L[oi]/Lfir = 1.6 x 10"^ with a standard 
deviation of about 50%, and the total range of values of roughly a factor of six (Figure 4) 
[O I] (145 iim) is fainter and was observed/detected only in a few galaxies. 




0.4 0.6 0.8 1 1.2 1.4 -1 -0.5 0.5 1 1.5 

F,(60 Mm)/F^(100 /xm) LoglLFm/LBiue) 



Fig. 4. — L[oi]/Lfir shows no trends with FIR colors or FIR-to-blue ratio, except for the 
non-detections for galaxies with [C II] deficiency. In detected sources alone, the observed 
L[oi]/Lfir ranges over a factor of six, as opposed to a factor of 60 for L[cii]/Lfir. 

The [O l](63/im) line is the other major coolant of neutral regions. Since the energy 
to excite this line is higher than for [C II] we expect this line to become more important 
relative to [C II] in warmer gas (i.e. higher Gq). This is indeed seen in Figure 5, where 
L[oi]/L[cii] increases with F,^(60 /im)/F,^(100 /xm), indicating that warmer gas correlates 
with warm dust. In Figure 1 we saw that L[cii]/Lfir decreases with F,^(60 //m)/F^(100 /im) 
and in Figure 4 we see that L[oi]/Lfir does not decrease with F^(60 //m)/F^(100 //m), so 
we should be able to predict the correlation seen in Figure 5. The generalized Kendall's 
rank test for the correlation between L[oi]/L[cii] and F,^(60 /im)/F,y(100 fxm) gives a 5.5cr 
significance for a sample of 50 galaxies. This is the strongest correlation we see in this 
dataset, and there is a good physical reason for that. PDR models (e.g. Kaufman et al. 
1999) show that both the increase in dust temperatures (F,^(60 iJ,m)/F^{100 nm)) and 
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[O I] (63 //m)/[C II](158//m) are due to an increase in FUV flux Gq. We make a detailed 
comparison between the PDR models and observations in the next section. 

[O I] (63 jjim) does not dominate the cooling for the galaxies in the sample 
where F,^(60 //m)/F,^(100 jim) < 0.8. Therefore, if photoelectric heating is the 
dominant heating mechanism, L[cii]/Lfir tracks the heating efficiency for the cooler 
galaxies [F,^(60 //m)/Fi,(100 /xm) < 0.8]. Total heating efficiency is given by the ratio 
(L[cii] + L[oi])/Lfir and decreases with Fi,(60 //m)/Fj,(100 //m), as discussed in the next 
section. 



2 



~ 0.8 
< 0.6 

0.4 
0.2 



0.4 0.6 0.8 1 1.2 

F^(60 /xm)/F^(100 ixm) 

Fig. 5. — The ratio of the main cooling lines for neutral gas, L[oi]/L[cii], shows a tight 
correlation with FIR colors (5.5o" significant), indicating that warm gas (which emits 
more [O I] (63/im) relative to [C II]) correlates with warm dust as indicated by high 
Fj,(60 /im)/Fj,(100 /im). The observed trend is consistent with the PDR models which 
ascribe the increase in temperature in gas and dust to higher FUV fiux Gq. [C II] dominates 
the coohng for galaxies with F,^(60 /im)/Fi.(100 /xm) < 0.8, and in galaxies with warmer 
dust [O l](63/im) becomes more important. 

The [O I] (63 jim) line is optically thick in some cases. This can be seen by 
comparing it to the [O I] (145 ^m) line. The maximum allowed ratio in the optically 
thin limit L[oi]i45/L[oi]63= 0.1 for T > 300K (TH85), whereas in the current sample 
L[oi]i45/L[oi]63=0.18 for NGC 3620. 

The L[oi]i45/L[oi]63 line ratio can also be used to probe gas temperature; L[oi]i45/L[oi]63 
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increases with gas temperature, but the S/N in the [O I] (145 fxm) hne is not enough in 
most galaxies where it is observed. 



4.4. The [O III] (SS/xm and 52/xm) lines 




0.4 0.6 0.8 1 1.2 1.4 -0.5 0.5 1 1.5 

F,(60 Mm)/F^(100 /j.ui) Log(L„B/LB,„J 



Fig. 6. — L[oiii]/Lfir does not vary with FIR colors (the correlation between L[oiii]/Lfir 
and F,^(60 fim)/F^{100 /im) is only 1.2 a significant). L[oiii]/Lfir shows an anticorrelation 
with Lfir/Lb (3.8 <j significant), which is largely contributed by the two irregular galaxies, 
NGC 1569 and IC 4662, which have the lowest Lfir/Lb and have very high Ljonij/LpiR- In 
fact, [O III] is the strongest far- infrared line in these galaxies. 

The [O III] line at 88/xm is one of the more easily detected lines. We detect this line 
in 24 out of 34 galaxies where it was observed. Two low luminosity, irregular galaxies, 
NGC 1569 and IC 4662, with low values of Lfir/Lb, have high values of L[oiii]/Lfir, with 
L[oiii] approaching or exceeding 1% of Lfir (Figure 6). The [O III] (88 /xm) line in these 
galaxies is stronger than the [C II] (158 /xm) and [O I] (63 /im) lines from neutral gas. There 
is other evidence that the line and continuum emission in these galaxies may be dominated 
by H II regions (see Hunter et al. 2001 for details). There is relatively little dust near these 
H II regions, hence the low extinction, low Lfir/Lb and high L[oiii]/Lfir- The irregular 
galaxies in the sample which are observed in [O III] have a high L[oiii]/Lfir- 

The [O 111] (52/xm) line is not so widely detected, or observed, due to the relatively 
low sensitivity of the detector at those wavelengths. This hne was detected in 3 out of 11 
galaxies where it was observed. 
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5. Interpretation of the data 

5.1. Electron densities derived from [O III] lines 

For three galaxies, NGC 5713, NGC 1569, and NGC 4490, we use the measured 
ratio of the two [O III] hues, at 88 and 52 //m, to derive electron densities Ue using the 
semi-empirical treatment of Rubin et al. (1994) with an assumed H II region electron 
temperature Tg = 10^ K. The results are fairly insensitive to the temperature assumed. 
The 52 fim upper limits in five other galaxies can be combined with the measurement of 88 
/im line flux to yield upper limits on ng (Table 3). 



Table 3: Ue from [O III] line ratios 



Galaxy Name [0 III 


88/ [0 III 


52 Log(ne/cm ^) 


NGC 0278 


> 1.75 


< 0.1 


NGC 520 


> 0.63 


< 2.6 


UGC 02855 


> 0.54 


< 2.7 


NGC 1482 


> 1.38 


< 1.8 


NGC 1569 


1.55 


1.5 


NGC 4418 


> 0.21 


< 3.4 


NGC 4490 


1.25 


2.0 


NGC 5713 


1.04 


2.2 



5.2. Variations in Ljcuj/LpiR and proposed explanations 

In Paper 1 we discussed many possible reasons for the extremely low values of 
L[cii]/Lfir observed in some galaxies. Some of these hypotheses can be ruled out with the 
observations of other FIR lines, the most prominent of these being [O I] (63 /xm) and the 
[Nil] (122 //.m). We review the various reasons proposed for the [C II] deficiency in Paper I 
and in Luhman et al. (1998) and discuss them in light of current observations. 

(A) The [C II] line could be optically thick in emission. This hypothesis is supported 
by the compactness of the mid-IR emission in the three galaxies with the lowest L[cii]/Lfir. 
The [C II] line becomes optically thick at A'"(C+) = 5 x lO^^cm"^ for a velocity width 
of 4 km/s (Russell et al. 1980). Assuming a velocity width of 120 km/s for NGC 4418 
(Sanders, Scoville & Soifer 1991) and N{C+)/N{H) ~ 1.4 x 10-^ optical depth of one in 
the [C II] line is reached for column densities of N[H) = 8 x 10^^ cm~^, or Ay = 40 in the 
PDR gas. The total extinction towards the nucleus of NGC 4418 is Ay > 50. This would 
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imply an unusually high ratio of PDR to cold molecular gas. 

The biggest uncertainty in estimating the column density of gas needed to produce 
optically thick emission is the velocity structure of the gas; if there is velocity crowding (e.g. 
near the center or in a bar) or if most of the self- absorption is local (i.e. in the same cloud 
or within a small velocity interval) , even a small column density of gas can be effective. But 
large optical depths arc difficult to produce because of the following argument. Models of 
PDRs, assuming that grain opacity and gas abundances both scale with metallicity, show 
that the column of in a PDR is only of order a few times lO"*^^ so that the [C II] has 
optical depths in a given cloud only of order unity, unless every line of sight through the 
cloud goes through many PDRs. But in that case one would see emission from the first 
optical depth of unity, so for an optical depth of 100 in C+, one would get 1 /100th of the 
emission compared to the optically thin case (and not e~^°°). So to explain L[cii]/Lfir 
which is 100 times lower, one would need an optical depth of 100 in the PDRs, which 
corresponds to a hydrogen column density of ~ 3 x lO^^cm"^. 

Another way to produce large optical depth is to have cold C"^ on the outside of the 
cloud while the PDR inside produces [C II] emission. The cold gas must then be turbulent 
enough to have a velocity width equal to or larger than the warm/emitting gas so that [C II] 
emission in the line wings does not escape. Self-absorption is seen in the 158/xm line when 
observed at high spectral resolution (Boreiko & Betz 1997), and in ISO observations of 
Sgr B2 (Cox et al. 1999), but the line width of the cold gas is usually smaller than that of 
the warm gas, leading to a reduction in fiux by a factor of roughly two, but cannot produce 
a factor of 60 reduction in flux without a good match in velocity widths. 

Since [O I] is present up to higher optical depths in clouds than C^, [O I] is more 
likely to have larger columns and to be self absorbed. Towards the Galactic center a deeper 
absorption proflle is seen in the [O I] (63 /im) line (Keene et al. 1999, Cox ct al. 1999, 
Poglitsch et al. 1996, Kraemer, Jackson & Lane, 1996) than in the [C II] line (Cox et al. 
1999). [O I] should go optically thick /aster than [C II] if optical depth were the reason for 
decreasing L[cii]/Lfir. Thus we should have seen a more dramatic decline in L[oi]/Lfir 
than is seen in L[cii]/Lfir with increasing F,^(60 //m)/F,^(100 //m) or Lfir/Lb. Since no 
decline is seen in L[oi]/Lfir (Figure 4b), we discount self-absorbtion as the dominant reason 
for the decrease in L[cii]/Lfir. 

(B) High dust extinction: For dust extinction to be important Tdust = 1 at 158 /im 
imphes Ay — 15, 000. This is still insufficient to explain the [C II] deficiency and is unlikely 

even for very obscured starbursts like Arp 220 (Fischer et al. 1999). Since the [O I] line 
at 63/im is at shorter wavelengths it is even more vulnerable to extinction effects but still 
requires far too high an Ay for this to be a viable explanation. We should have also seen 
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a dramatic decrease in L[oi]/Lfir with increasing Lfir/Lb in figure 4b since Lfir/Lb is a 
measure of extinction effects in galaxies, but such a trend is not seen. 

(C) Softer radiation from older/less massive stellar populations could lead to lower 
L[cii]/Lfir due to a lack of C-ionizing photons. If much of the grain heating is by longer 
wavelength photons, rather than by FUV, [C II] will be relatively weak compared to FIR 
because grain photoelectric heating of the gas will be inefficient. This can happen when 
dust heating is by old, low mass stars. Luhman et al. (1998) postulate this scenario for 
[C Il]-deficient ULIRGs. An aging starburst and a weak UV field relative to a softer, 
optical/infrared field which heats the dust, is also used to explain the [C II] line weakness 
in the Galactic center (Nakagawa et al. 1995). This mechanism seems unlikely to explain 
the most dramatic [C II] deficiency systematically found in galaxies with warmer FIR colors 
and higher Lfir/Lb ratios, with the hotter dust pointing to the presence of massive stars 
and UV photons. It also seems unhkely that all starbursts are old by the time we see them. 
This explanation is more likely for the decrease seen in L[cii]/Lfir in early type galaxies 
with low rates of star formation and the lowest Lfir/Lb in the sample (to the left extreme 
of Figure lb; also see Malhotra et al. 2000). Softer radiation fields in early type galaxies 
then lead to somewhat lower (by factors of a few) L[cii]/Lfir; this is corroborated by lower 
values of [C II] /CO seen in one of these galaxies (Malhotra et al. 2000). 

(D) An obscured AGN can lead to both higher Lfir/Lb and lower L[cii]/Lfir- This 
hypothesis is supported by the observed compactness of the [C Il]-deficient sources in the 
mid-IR (Dale et al. 2000). A buried AGN could produce a small L[cii]/Lfir since the UV 
field is inefficient at making [C II], both due to its hardness (higher ionization states for C 
will be common) and due to its overall strength (see (F) below). This explanation may not 
hold for the ULIRGs, e.g. Arp 220, which are now believed to be powered by starbursts 
(Genzel et al. 1998). The AGN hypothesis may explain the deficiency of [C II] but does not 
easily explain the trends in L[cii]/Lfir with Lfir/Lb and F,^(60 //m)/F,^(100 /im) unless a 
significant fraction of galaxies have obscured active nuclei, and even then one would expect 
the AGNs to add more scatter than a trend. 

(E) For high values of Gq (or T) and n, [O I] (63//m) emission is stronger than [C II] 
emission, since [C II] is quenched by coUisional de-excitation. The grain charge, which 
determines the gas heating efficiency as measured by the observed ratio (Lpn] + L[oi])/Lfir, 
is set by Go/n. Therefore, with Go/n fixed, raising Gq and n together does not change 
(Lfcii] + L[oi])/Lfir, but it does increase L[oi]/L[cii] and L[oi]/Lfir while decreasing 
L[cii]/Lfir. In other words, the heating efficiency stays constant, but the cooling 
increasingly emerges in the [O I] line instead of the [C II] line as Go rises with Go/n 
fixed. In section 4.3, we observed that the ratio of the two major coohng fines of the 
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neutral ISM, L[oi]/L[cii], increases with increasing F,^(60 fj,m)/F^{100 //m) (Figure 5). 
F,^(60 ;um)/Fj,(100 /xm) is a rough indicator of the dust temperature and increases with the 
FUV flux Go- Since the [O I] (63/im) hue has an excitation energy AE//i;=228 K compared 
with AE//c=92 K for the [C II] (158/im) hne, an increase in gas temperature caused by an 
increase in Gq will increase Ljoij/Lpiij. Therefore, the correlation between L[oi]/L[cii] and 
Fi,(60 /i.m)/Fi,(100 /im) is likely due to an increase in Gq, which increases both gas and dust 
temperatures. However, Figure 7 shows that (Ljcn] + L[oi])/Lfir does not stay constant 
with increasing Fi,(60 /xm)/Fi.(100 //m) or increasing Gq. This is a clear indication that the 
density is not rising as rapidly as Gq and that grains are becoming more positively charged. 
Therefore this hypothesis, although it may contribute somewhat to the [C II] deficiency, 
cannot be the dominant cause. 

(F) For high ratios of FUV flux to gas density (Go/n), the dust grains become highly 
positively charged and are less efficient at heating the gas because of a higher potential 
barrier to photoelectric ejection (TH85). Hence there is less total gas cooling relative to the 
FIR continuum emitted by grains. For example, increasing Gq while keeping n fixed raises 
the grain charge, lowers the gas heating efficiency (L[cii] + L[oi])/Lfir, lowers L[cii]/Lfir 
and raises L[oi]/L[cii]. An increase in Gq, accompanied by an increase in Go/n, is the likely 
explanation for the trend of decreasing (Ljcn] + L[oi])/Lfir with F,^(60 /xm)/F,^(100 /xm) 
seen in Figure 7, as the dust gets warmer with increasing Go and more positively charged 
with increasing Gq/u. In addition, it accounts for the decrease in L[cii]/Ltir with 
Fi,(60 //m)/Fi,(100 //m) seen in the data plotted in Figure 9. This, then, is our favored 
hypothesis; in the next section we compare PDR models to L[cii]/Lfir to further bolster 
this hypothesis. 



5.2.1. [Cll] deficiency Case Study: Arp 220 

Of the explanations suggested above, (C) & (F) postulate a decrease in neutral gas 
heating relative to grain heating. In the others, either [C II] is attenuated by dust (B) or 
the coohng emerges in other fines (A, D and E). We have already ruled out E (cooling by 
[O I]) as the dominant cause of the reduction in L[cii]/Lfir. We can further test hypotheses 
A and D by measuring the emission in all the major cooling lines to see which, if any, 
line has replaced [C II] and [O I] as the coolant of the neutral medium. The total cooling 
budget can then be computed to see if the heating efficiency has declined. (The test is 
less stringent for D, since the gas cooling may be predominantly in optical lines which are 
absorbed by dust and emerge as continuum.) 

At present Arp 220 is the only [C II] deficient galaxy for which a relatively complete 
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Fig. 7. — Heating efficiency (L[cii]+L[oi])/Lfir vs grain temperature Fi,(60 //m)/Fj,(100 fxm). 
A decline in the total cooling in these lines with respect to FIR is seen in galaxies with 
warmer dust temperatures as indicated by Fi^(60 ;um)/Fi^(100 /xm). This is expected as 
grains become more positively charged in regimes of high Gq/h. The rank correlation test 
yields that the anticorrelation between (L[cii] + L[oi])/Lfir and Fi,(60 /xm)/F,^(100 /xm) is 
3.3(7 significant. 

census of gas coohng lines has been made. Gerin and PhiUips (1998) have observed this 
galaxy in the [CI] (492 GHz) line and higher-lying CO rotational lines, and Sturm et al. 
(1996) observed this galaxy in H2 rotational lines. Table 4 summarizes the energy output 
in cooling lines of neutral gas in Arp 220. 

We see that only the molecular hydrogen lines come close to matching even the reduced 
output in [C II]. [O I] (63/xm) is seen in absorption (Fischer et al. 1999). We conclude that 
there is a reduced ratio of gas heating to grain heating in at least this [C II] deficient galaxy. 
This is then consistent with hypotheses C and F. 



5.3. Trends in L[cii]/Lfir: comparison to the PDR models 

Of the various explanations for the trends and the deficiency of L[cii]/Lfir listed above, 
the data most favor scenario (F) where the decrease in L[cii]/Lfir is due to reduced heating 
via charged grains at higher values of Gq/u, with scenario (E) contributing to the decrease 



- 21 - 



Table 4: Arp 220: Energy Budget 



Flux (observed) [W/m ^) Flux(extinction corrected) {W/m ^) 



FIR Continuum ^ 


6.8 X 10^12 




[C II] 1 


8.7 X 10-16 




CI i'Pi-'Po) ' 


1.9 X 10-1^ 




CI (total) 1 


1.0 X 10-16 




CO (1-0) 1 


1.8 X 10-1^ 




CO (total) 1 


2.1 X 10-16 




i/2 S(5) 2 


2.4 X 10-16 


2.4 X 10-15 


i/2 S(2) 2 


< 1.5 X 10-16 


< 7.0 X 10-16 


S(l) 2 


2.3 X 10-16 


9.7 X 10-16 


i/2 S(0) 2 


< 3.5 X 10-16 


< 7.3 X 10-16 



1 Gerin & Phillips 1998 

2 Sturm et al. 1996 

in [C II] luminosity when [O l](63 fim.) starts dominating (F,^(60 ;um)/F,^(100 fira) > 0.8). 
In Figure 9a and 9b we see the quantitative comparison of L[cii]/Lfir and L[oi]/Lfir 
with the PDR models of Kaufman et al. (1999). The term PDR refers to the neutral 
ISM in galaxies, including diffuse clouds or cold neutral medium (CNM), warm neutral 
medium (WNM), dense neutral gas around H II regions and molecular clouds illuminated 
by the interstellar radiation field (ISRF). PDRs are formed wherever FUV photons play a 
dominant role in the heating and chemistry of the gas. 

In order to compare the data to the models we need to correct the quantity Lpm 
used previously to include the dust emission outside the wavelength range 40-120/xm. 
Since the long wavelength part of the dust emission (A > 120//m) is not measured for 
these galaxies, we have to rely on dust emission models based on Galactic dust and other 
well-studied galaxies. Empirical data from IRAS and ISO indicate that the overall spectral 
shape in the infrared may be represented to first order as a function of one parameter, the 
F,^(60 /im)/F^(100 fim) ratio. Dale et al. (2001) have proposed a single parameter family of 
spectral energy distributions (SEDs) based on those data from which we convert FIR (dust 
emission in the 40-120/xm) to total infrared emission (TIR; 3-1100 /xm), including emission 
from PAHs and dust associated with cold cirrus. The uncertainty in converting FIR to TIR 
is about 15% (D. Dale, private communication). 

Two things are apparent from this comparison. First, the PDR models do reproduce 
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Fig. 8. — Comparison of L[cii]/Ltir with PDR models. TIR is the total IR emission (3-1100 
^m; see text). The solid dots represent the L[cii]/Ltir values for the 60 galaxies in this 
sample. The curves represent L[cii]/Ltir values for gas densities from 10^ to lO^cm"^, with 
notches indicated for the progression of Go, from the PDR models of Kaufman et al. 1999. 
Lines of constant Go are vertical in this plot. The FIR colors become warmer for higher Go- 
We see that the trend of decreasing L[cii]/Ltir with Fi,(60 /xm)/F,^(100 /xm) is explained by 
the models as being due to the increase in Go. 

the decrease in L[cii]/Ltir with Fj,(60 /im)/Fj,(100 //m). For the extremely deficient 
L[cii]/Ltir galaxies, the PDR models of Kaufman et al. (1999) require Gq/u > 10^ cm^ 
in order for L[cii]/Ltir < 10""^. We note in Figure 8 that a second PDR solution for 
L[cii]/Ltir < 10~^ occurs for a much higher Go/n ratio. However, this solution produces far 
more L[oi] than is observed; it is essentially hypothesis (E) ruled out from the last section. 
In the regime Gq/u > 10^ cm'^, radiation pressure on grains may drive them through the 
gas at large drift velocities, resulting in additional gas heating, a detail not considered in 
the Kaufman et al. (1999) models; recent calculations by Weingartner & Draine (1999) 
show that such grain drift will only be important for large grains, while small grains, which 
are responsible for gas heating and much of the FUV absorption in the PDR models, will 
not be significantly affected. 

Second, we note in Figure 8 that about half the galaxies have L[cii]/Ltir which is 
too high by a factor of about two to be explained by PDR models with any parameters. 
This is not very surprising, since not all of the [C II] comes from PDRs; ionized gas in 
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diffuse ionized regions can contribute significantly to [C II] emission as well (Petuchowski & 
Bennett 1993, Heiles 1994). We discuss this further in the next section. 



5.4. [C II]/[NII](122 /im) from ionized regions 

The ionized gas contribution to [C II] can be estimated from measurements of the 
[Nil] (122 fjLm) line, which arises exclusively in the ionized gas. For the asymptotic high 
electron density limit {ue » ricrit, where ricrit = 3.1 X lO^cm-^ for [Nil] (122 /xm) and 
ricrit = 50cm-3 for [C Il](158//m)) (see Table 1): 

^l-l- 0.288-'^(^"' 



L|Niil N(N+) 

where N{N~^) and N{C~^) are the numbers of the respective ions, and the Einstein 
A coefficients used are Ai^s = 2.3 x 10~^s~^ (Froese-Fischer & Saha 1985) and 
Ai22 = 7.4 X lO-^s"^ (Froese-Fischer 1983). 

In the low density limit (rig << ricrit) the ratio is given by: 

L[cii] _. N{C^) 
L[Nii] ■ N{N+) 

where the collision strengths are taken to be = 2.15 (Blum and Pradhan 1992) and 
1^13 = 0.272 (Lennon and Burke 1994). 

The ratio of the [C Il](158/xm) and [NII](122/xm) luminosities depends on the gas phase 
abundances of these ions. We expect the and A^"*" abundances to be different in different 
galaxies in the sample. At present, the abundances in these galaxies are not known so we 
use Galactic (Solar, or interstellar) abundances. Even in the Galaxy the differences are 
substantial: [C/H] = lA x 10~^ in absorption line measurements of diffuse gas (Sofia et 
al. 1997) and [C/H] = 2.4 x 10"^ (Esteban et al. 1998, Rubin et al. 1993) as measured 
in emission in Orion. Nitrogen abundances are less discrepant: [N/H] — 7.5 x 10~^ from 
absorption fine measurements of the diffuse gas (Meyer et al. 1997) and [N/H] — 6.8 x 10~^ 
from emission line studies in Orion (Esteban et al. 1998, Rubin et al. 1993). For various 
combinations of those abundances, the expected L[cii]/L[nii] ranges from 5.7 to 10.7 in 
the diffuse {n « ricrit) limit and from 0.54 to 1.0 for dense ionized gas (n >> ricrit)- 
Furthermore, combining the abundances observed in the diffuse ionized gas (DIG) with the 
diffuse gas estimates points to 5.7 as the most reasonable value to adopt. 

In the galaxies where both lines are detected, L[cn]/L[Nii] ranges from 4.3-24, and 
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including upper limits increases the range to 4.3-29. The geometric mean of L[cii]/Lnii ratio 
in the sub-sample where both lines are detected is 8, and the arithmetic mean in 9. 

The total [C II] in a galaxy comes from PDRs, DIG and H II regions, whereas the 
[Nil] arises from DIG and H II regions. Thus the variable fractions of the galaxy's ISM 
involved in each of those components will add scatter to Ljchj/Ljnh], as will variations in C 
and N abundances. Using a constant multiplied by the [Nil] flux to correct the [C II] flux 
for ionized region contributions is therefore an oversimplification we adopt in the absence 
of further information on individual galaxies. The observed values of L[cii]/L[nii] cluster 
clearly near the low end, with more than half between 4.5 and 6, suggesting that the most 
generally applicable correction to [C II] must be near the lower end of that histogram as 
well. This also favors the ratio of 5.7 based on the lower abundance of C+, which we adopt 
in what follows. 

The amount of [Nil] emission from diffuse versus dense regions has only been measured 
for the Milky Way. From FIRAS measurements of the two [Nil] lines at 205 and 122/xm in 
the Milky Way, Bennett et al. (1994) conclude that most of the emission is from diffuse gas. 
The predicted ratio of these lines is I(122)/I(205)=0.7 for diffuse gas {ue « 100 cm~^), 
and I(122)/I(205)=3 for Ue ~ 100cm-=^ (Rubin 1985). Taking the middle of the range 
of the observed ratio I(122)/I(205)=1.0 to 1.6 (Wright et al. 1991), about 75% of [Nil] 
emission comes from diffuse gas. This implies that [C II] emission from the ionized medium 
[CII]ion = 4.3 X [iV//](122). We will therefore approximate the [C II] flux from PDRs as 
[C//]c = [CI I] - 4.3 X [iV//](122) in what follows. 

The mean and median adjustment made in the [C II] flux to thus remove the ionized gas 
contribution is roughly 50%. This estimate agrees with a similar estimate by Petuchowski 
& Bennett (1993) for the contribution to the Milky Way [C II] by the warm ionized 
medium. This is not surprising since we have used Galactic abundance ratios, and the 
average L[cii]/L[nii] hi this extragalactic sample agrees with the Milky Way value. We 
believe that while statistically this may be a reasonable treatment to estimate the ionized 
gas contribution to [G II], it may be inaccurate for at least some individual galaxies because 
of variations in G/N abundance ratios and in the fraction of diffuse and dense ionized gas. 
In Section 6.1, we check for and discuss the accuracy of estimating [G II] c. 

5.5. [O I]/[C II] compared to FDR models 

In section 4.3, we noted a strong correlation between the ratio of the two PDR cooling 
lines [O l](63 //m)/[G Il](158 //m) and far-infrared colors, F^(60 //m)/F^(100 //m). The 



-25- 




0.4 0.6 0.8 1 1.2 ' 0.4 0.6 0.8 1 1.2 

F^(60 jU.m)/F,(100 /xm) F^(60 Atm)/F„(100 fu.m) 



Fig. 9. — (a) The observed correlation between the ratio of the two PDR coohng hnes 
L[oi]/L[cii] and the dust temperature as indicated by Fi,(60 /im)/Fi,(100 iim) is shown by 
the sohd dots and upper hmits. This is compared to theoretical PDR models (Kaufman et al. 
1999). While the models run parallel to the trends seen, L[oi]/L[cii] in some models is too high 
for about half the sample, (b) The same comparison is made, except we attempt to subtract 
the contribution to [C II] luminosity from ionized gas by scaling with [Nil] luminosity. The 
solid points are where [Nil] was detected, and the triangles are where we used the 2-a 
upper limits on the [Nil] measurements. We see that the scatter in the correlation increases 
because of the uncertainties in [Nil] observations and the scaling factor used to subtract the 
[C II] emission due to ionized gas (cf section 5.1). The correlation between L[oi]/L[cii] and 
Fj,(60 yLim)/Fi,(100 /xm) persists, and the PDR models better fit the data. 

increase in L[oi]/L[cii] is due to the drop in L[cii] compared with LpiR and not due to an 
increase in L[oi]/Lfir. 

Qualitatively, this trend can be explained by postulating that both the increase in 
L[oi]/L[cii] and Fi^(60 /im)/F,^(100 /xm) come about in PDRs as a result of increase in 
the average FUV flux Gq. The overall decrease of (L[cii] + L[oi])/Lfir (Figure 7) means 
that Go/n is also increasing. In PDR models, increasing Gq leads to an increase in the 
grain temperature and therefore an increase in Fj,(60 iim)/F^{100 /im). Increasing Go also 
raises the gas temperature which raises the ratio L[oi]/L[cii]. Figure 9 shows a quantitative 
comparison between the observations and the PDR models of Kaufman et al. (1999). While 
the models run parallel to the trends seen, the observed L[oi]/L[cii] is too high compared to 
the models. Some of this could be due to [C II] flux from ionized gas. 
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When the contribution from ionized gas is subtracted using [Nil] (122 //m) line 
measurements as discussed in the last section, the models and the data show better 
agreement. In addition, by using the corrected, [C II]c, in L[oi]/L[cii], we see that the 
increase in L[oi]/L[cii] is due to an increase in both FUV flux. Comparison with PDR 
models (Figure 9b) shows that for n = 10 — lO^cm^^ the increase in L[oi]/L[cii] is due 
to an increase in Gq, but for n > lO^cm^^ the L[oi]/L[cii] increases due to an increase 
in both Go and n. The critical density for [C II] is 3 x lO^cm"^, while for [O I] it is 
5 X 10^ (T/300) cm~^. It is interesting, then, that L[oi]/L[cii] shows a good correlation with 
the FIR color, which does not depend on n but on Gq alone. In the next section we will 
present more evidence that Go and n tend to increase together. 

6. Physical Conditions in the PDRs 

After subtracting an estimated fraction of [C II] line flux that arises in ionized gas, we 
are left with [C II] line flux arising from PDRs. To derive the physical conditions of PDRs we 
compare the observed line ratios and line to continuum ratios with a grid of recent detailed 
PDR models by Kaufman et al. (1999) which take into account the chemistry, heating and 
cooling in PDRs. These models calculate the line emission in [C Il](158//m), [O I] (63 //m 
and 145 /^m), and dust continuum emission for a plane slab of gas illuminated from one 
side by FUV. Gas heating is dominated by photoelectrons ejected from classical (i.e. big, 
spherical) grains, and from PAHs following the treatment by Bakes and Tielens (1994). 
Including the effects of PAH photoelectric emission raises the PDR surface temperatures by 
as much as a factor of 3 compared to the previous PDR models (e.g. WTH). 

The comparison of measured fluxes of [C II] (158/xm) and [O I] (63 fim) lines to the 
PDR model predictions allows us to infer the physical conditions in the PDRs, primarily 
gas density n and FUV flux Go. This method may not give the right answers if: (1) the 
various lines and the continuum arise from different locations in the galaxies; (2) if there 
is self-absorption of the lines from cold gas outside the PDRs where they originate. [O I] 
(63 /im) is most susceptible to this, and we try to correct for it in the models as described 
below. 

The PDR models are used to calculate the emergent flux from the front face of a plane 
parallel slab of gas illuminated from one side. A galaxy has many PDRs at all orientations, 
and optical depth effects are non- negligible. In the approximation that most dense PDRs 
are the shells of molecular clouds, and that the [C II] line and dust continuum emission 
are optically thin whereas [O I] (63 /xm) is optically thick, [O I] is seen only from the 
front side of each cloud and [C II] and FIR from both the front and the back sides. The 
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velocity dispersion from cloud to cloud however allows most [O I] (63//m) photons that 
have escaped their parent cloud to escape the galaxy entirely. This scenario implies that 
we should observe only half of the [O I] flux and all of the [C II] and dust continuum flux 
expected from PDR models. This adjustment has been applied hereafter to the PDR model 
predictions. 

Figure 10(a) shows the measured ratios Ljoij/Ljcn],, and (Ljoij + LpujJ/LxiR compared 
with the values obtained from a grid of PDR models with various values of Gq and n. With 
just these two measured ratios there are two possible regimes of Go and n which reproduce 
the observed values. One of the regimes is shown in Figure 10(a), while the other is a high 
n (10*^cm~^ <n < lO^cm"^), low Gq (Gq — 1) regime. To distinguish between the two, we 
need to reliably measure the ratio L[oi]i45/L[oi]63- This measurement is not available for 
most of the galaxies in this sample, [O I] (145 /xm) being too faint. However the high n, low 
Go solution does not reproduce the line fluxes measured, i.e. the fluxes in the [C II] line 
fall short by about 2 orders of magnitude when 10^ < n < 10^ and Gq ~ 1 compared to 
the observed values. This is not to say that such high density, low radiation environment 
(e.g. molecular clouds) do not exist in galaxies, they simply do not contribute much to the 
observed fluxes seen in [O I] , [C II] lines and the FIR continuum. 

6.1. Go and n from [O l]/TIR values 

We can also use the L[oi]/Ltir ratio to derive Gq and n by comparing the observed 
L[oi]/ Ltir vs. F,^(60 //.m)/F^(100 fim) with models. Such a comparison, shown in 
Figure 11(b), using models by Kaufman et al. (1999), suggests that the constancy of 
[O I]/F1R (or TIR) against variations in F,^(60 ;um)/F,^(100 /im) implies a simultaneous 
increase in both gas density n and FUV flux Go- The model Fi,(60 /xm)/Fi,(100 /xm) ratio 
for a given UV flux Gq is calculated following the treatment by HoUenbach, Takahashi 
& Tielens (1991). The 60 and 100 //m fluxes are calculated assuming classical grains in 
thermal equilibrium. For low Go, the emission from classical grains is from the Wein side 
of the blackbody curve, and the Fi,(60 //m)/Fi,(100 /xm) is set by stochastically heated 
grains. This has the effect of yielding essentially the same F;^(60 /im)/Fiy(100 ^m) for 
Go = 1 — 10^'^. The other uncertainty in using L^qi]/ Ltir to derive Gq and n is that [O I] (63 
/im) may be optically thick; we have accounted for this by reducing L^qi]/ Ltir from the 
models by a factor of two, as before. An advantage to this disagnostic however is that we 
avoid uncertainties associated with correcting for the [C II] flux from ionized regions. 

Figure 11 compares the Gq and n values obtained using each set of observables, 
L[oii/L[cii]c vs (L[oi] + L[ciiiJ/Ltir and L^oi]/Ltir vs F^(60 /xm)/Fj,(100 /xm). The 
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Fig. 10. — (a) The measured ratios L[oi]/L[cii]c (L[oi]+L[cii]c)/Ltir (circles and triangles) 
are overplotted with a grid of PDR model values of the same. Solid lines are contours of 
constant gas density n and the dotted lines are contours of constant FUV radiation flux 
Go- Two regimes of Gq and n reproduce the observed ratios. Only one is shown here 
for clarity. The other regime is ruled out as it would not produce the observed [C II] 
and [O I] fluxes. The observed points are as in Figure 10: filled circles represent values 
where [C II] was corrected for the contribution from ionized gas using measured values of 
[Nil], and triangles are where upper limits for [Nil] were used, (b) The measured ratios 
of L[oi]l Ltir and Fy(60 /im)/F,^(100 /im) are plotted along with a grid of L[oi]/Ltir and 
Fi,(60 fim) /Fi,(100 /xm) values calculated according to PDR models of Kaufman et al. (1999) 
and dust models of HoUenbach, Takahashi & Tielens (1991). 



difference in log(Go) obtained by these two methods ALog{Go) shows a mean of 0.08 
and Alog{n) shows a mean of 0.06. Thus there is no indication that one method yields 
systematically higher or lower values of Go or n. From the absolute values of AGo and 
An, we infer that for half the galaxies Go and n values from the two methods agree to 
better than a factor of 2 and 1.5, respectively. This difference may be taken as an estimate 
of uncertainty in the derived values, which is sensitive to some measurement errors (one 
method uses only the [O I] line while the other uses [C II], [O I] and [Nil]) and to some 
sources of systematic errors (for example, subtraction of the ionized gas contribution to 
[C II] flux), but not to all possible systematics. For instance, this estimate would not 
be sensitive to uncertainties in the abundances used in the models, or to the assumption 
that about half the dust heating is due to FUV photons, or to the fact that we measure 
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luminosity weighted averages of line and continuum contributions emanating from many 
parcels of ISM with different physical conditions. 




Log(n) 

Fig. 11. — This figure shows the Go and n solution for galaxies based on comparison of 
data and PDR models shown in figure 11. The filled circles are Go and n values estimated 
from L[oi]/L[cii]o and (L[oi] + L[cii]c)/Ltir; and the open squares are the Go and n values 
derived from L[oi]/Ltir and F,^(60 /xm)/Fi^(100 /im). The Gq and n values derived from 
the two methods are connected for each source to give an estimate of the uncertainty in the 
parameters. A least squares fit is made to both sets of Gq and n values assuming equal error 
in both axes. The best fit slopes are 1.4 and 1.3 respectively, i.e. Gq oc n°, with a = 1.3 — 1.4, 
which is consistent with the emission coming from PDRs surrounding ionization-bounded 
expanding H II regions. 



6.2. Discussion and interpretation 

The average physical conditions in this sample of normal galaxies are FUV flux 
Go = 10^ ~ 10^'^ and gas densities n = 10^ — lO^'^cm"^. The derived values of Gq and 
n follow a trend with Go cx: n", a = 1.3 or 1.4 (Figure 12). This may be evidence for a 
Schmidt law on fairly local scales in the galaxies, where high gas densities correlate with 
high star formation rates (Schmidt 1959) and agrees with the recent determinations by 
Kennicutt (1998) where the surface densities of star-formation and gas scale by an exponent 
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N=1.4: T^sFR oc ^I'as- The observed scaling between Gq and n may be also dictated simply 
by geometry if the FIR line and continuum emission is dominated by regions near young 
stars. From simple Stromgren sphere calculations (cf. Spitzer 1978) we can derive that 
the FUV flux at the neutral surface just outside the Stromgren sphere should scale as 
Go oc n^^^, which is consistent (within errors) to the scaling seen in Figure 11. 

We favor the latter interpretation primarily because it is the more direct one. The 
Schmidt law applies to surface density or a density which is averaged over the entire galactic 
disk thickness. It reflects the physics that on a global scale, the star formation rate may 
scale as the mass (density) or, if cloud collisions are important, for example, on a higher 
power of the mass (density). However, the high values of Gq and n which we derive make 
it clear that we are probing a very local phenomenon. We are not probing the rate of star 
formation in as much as the typical distance from an OB star or OB association to a PDR 
which absorbs its FUV photons. O and early B stars form in dense cores within giant 
molecular clouds (GMCs), whose diameters are of order 10-50 pc. For a fraction of their 
lives, the OB stars are embedded in the clouds, surrounded by very dense and compact 
H II regions and PDRs illuminated by intense UV fluxes. Later, the expanding H II regions 
break through the surfaces of the clouds, champagne flows are created, and eventually the 
OB stars lie outside, but close, to their natal GMCs. To produce a flux of Go ~ 1000 on the 
surface of that cloud, an 07 star must be only about 1 pc distant. With a relative speed 
of 1 km/s with respect to the cloud, an O star travels only 1 pc in 10^ years, a substantial 
portion of its lifetime. Therefore, the FUV photons from OB stars tend to be absorbed in 
two very different environments. A substantial fraction (McKee & Williams 1997 estimate 
> 30% for the Lyman continuum photons) are absorbed in their immediate vicinity - in the 
GMC of their birth. The rest escape from the natal cloud and travel through the diffuse 
ISM to be absorbed primarily by diffuse gas at distances of order 100-300 pc, with much 
lower values of Go and n (~2 and 100 cm~^ in the local ISM of the Milky Way). These two 
components are consistent with the decomposition of the FIR dust emission into "active" 
and "cirrus" components by Helou (1986). [C II] (158/im) can have significant contributions 
from the DIG, the diffuse neutral gas, and the locally irradiated PDRs of the natal GMCs. 
[O I] (63/im) is dominated by emission from the PDRs of the natal GMCs. The PDR 
modeling procedure provides an average of these two components, weighted to the high 
Go and high n GMC component. The simple Stromgren sphere scaling Go oc n^^^ applies 
directly to this GMC component, and should therefore manifest itself in the integrated 
fluxes as seen in Figure 12. 

The derived PDR pressures Q also lend credence to the view that the Go oc n^'^ 



The PDR pressure is calculated as P = nkT where gas density n and temperature T are derived from 
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correlation arises from the correlation of n and Gq in an expanding H II region surrounded 
by a PDR. In such a picture, the H II region thermal pressure should approximately equal 
the PDR thermal pressure. The thermal pressure in the H II regions can be estimated 
from the electron density derived from [O III] (88;Um)/[0 III] (52yum) ratio if both lines 
are observed (Rubin et al. 1994). This comparison was made for the three galaxies in our 
sample with measurements of the [O III] lines. The pressures derived for H II regions are 
tabulated in Table 7. The H II region and PDR pressures agree to roughly a factor of two 
in all cases (assuming Tg = SOOOX). 
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Fig. 12. — The average pressure near star forming regions in the galaxies is plotted against 
the FIR colors. The pressure is a product of the gas densities n and temperatures T 
derived from the PDR models. Both n and T increase with Fi,(60 /im)/F,^(100 /xm), but the 
increase in n is more dramatic and is responsible for the correlation between pressure and 
Fi,(60 //m)/Fi,(100 jim). These pressures likely do not reflect the average thermal pressure 
of the ISM, but do represent the pressures in PDRs surrounding H II regions associated with 
O and B stars. Curves are shown for the pressure in PDRs surrounding H II regions formed 
by stars of different spectral types (see text). 

The derived temperatures at the PDR surfaces range from 270-900 K, and the pressures 
range from 6 x 10^ — 1.5 x lO^Kcm"^ (Table 6). The lower value of the pressure range 




PDR models (see Kaufman et al. 1999 for details) 
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is roughly twice the local solar neighborhood value and the upper end is comparable to 
pressures in H II regions in starburst galaxies (Heckman, Armus & Miley 1990) which also 
corresponds to the pressure and surface brightness at which starbursts saturate (Meurer et 
al. 1997). 

Figure 12 shows that the thermal gas pressures increase with grain temperatures, or 

the flux of photons incident upon the absorbing grains. This is mainly due to the L[oi]/L[cii] 
vs Fi^(60 ;Lim)/Fi^(100 /xm) correlation in Figure 5. Higher values of L[oi]/L[cii] imply higher 
Go and gas densities, and therefore higher pressures. The temperature of the gas also 
increases with F,^(60 /im)/F^(100 fim) and contributes to the increase in pressure, but the 
increase in the gas density n is the more dramatic of the two. This too is explicable in terms 
of the local absorption of photons by the expanding H II regions and associated PDRs in 
the natal GMCs of OB stars. Whether grains absorb most of the flux in the H II region or 
in the neighboring PDR, the photons will be absorbed at roughly Rs, the Stromgren radius 
for that star. If 0, is the Lyman continuum photon luminosity of the star, then the photon 
flux is proportional to (pi/Rs oc nl^^ oc Pyn = -PpdR) assuming pressure equilibrium between 
the H II gas and the PDR regions associated with them. Since the grain temperature scales 
roughly as the photon flux to the 1/5*^ power (HoUenbach, Takahashi & Tielens 1991), 

Igr CK. -TpDR- 

We show a simple model to help explain the observed correlation between the FUV 
field strength, Gq, and the PDR pressure, P. If we assume that PDR gas is illuminated 
by hot stars in star clusters, then there is an approximate natural scaling between the 
FUV fiux from stars in the cluster, the H II region density, and the pressure at the H II 
region/PDR interface. Assuming that all of the FUV flux emerges from one star of a given 
spectral type which has a bolometric luminosity, Lf,oi and which produces S — 10^^ 5*49 
photons s~^ of ionizing photons, the size of the surrounding H II region, Rs is set by the 
number of ionizing photons and the electron density in the H II region. Assuming further 
an electron temperature, T,. = 10^ K, the combination UeTe sets the H II region pressure. 
Then, for a given bolometric stellar flux, the value of Gq at the PDR surface is determined 
by Go = Li,oi/4:TtR'^. If there is a pressure equilibrium at the H II region/PDR interface, 
then we find a correlation between Go and P. This is plotted as a correlation between 
F;^(60 fj,m)/F^{100 fj,m) and P (Figure 13). The calculations use three different spectral 
types (05, 09, and B3) and electron densities ranging from 10^ to lO^cm"^. Zero-age main 
sequence values of L^oi and 5'49 for the various spectral types are from Parravano (private 
communication) . 

Because the ISO-LWS beam views an entire galaxy for our sample galaxies, it 
incorporates many tens to hundreds of thousands of O and early B stars, with their 
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associated H II regions and PDRs. As discussed above, the Gq, n and P we derive for 
a given galaxy represent an average value, weighted to the dense natal GMCs which lie 
close to the OB stars. The range of the Gq, n and P for the different galaxies can reflect 
several interesting differences in their star formation processes and histories. If global star 
formation in normal galaxies occurs in bursts, rather than continuously, then those galaxies 
which are observed shortly after their bursts few million years) will show high Gq, n and 
P because the OB stars will not have had time to travel very far from their natal clouds. 
This hypothesis, however, requires the star formation bursts to be of significant amplitude, 
resulting in sustantially different instantaneous mass distributions of ionizing stars among 
galaxies. 

The differences in Gq, n and P may instead reflect differences from galaxy to galaxy in 
the GMCs which form the OB stars. For example, larger GMCs may keep their OB stars 
embedded for a longer fraction of their lifetime, resulting in higher average Go and n. We 
note that n is the density in the PDR, which can be considerably higher than the average 
density in the GMC. On the other hand, the GMCs could be the same size but denser. In 
this scenario, the higher density ambient gas would lengthen the embedded phase of the 
OB star, and result in higher average Gq and n. In either case, the galaxies with higher 
derived Go and n would contain more massive GMCs, on average. 



7. Summciry and Conclusions 

In this paper we have attempted to understand the energetics and the physical 
conditions in a statistically representative set of star-forming normal galaxies by studying 
the atomic and ionic fine structure lines in the far-infrared. Such sensitive observations 
of a large sample were made possible by having a cryogenically-cooled observatory, ISO 
(Kessler et al. 1996), in space. The sample was selected to span a range in properties 
such as morphology, FIR colors Fi,(60 //m)/Fj,(100 /xm) (indicating dust temperatures), 
and Lfir/Lb (indicating star-formation activity and optical depth). For a randomly drawn 
sample of galaxies, many of these parameters are correlated. Care was taken to span the 
full range of parameter space as much as possible in this sample. Still, there remains some 
correlations between these parameters. Galaxies with high Lfir/Lb also tend to have 
warmer FIR colors and tend to be luminous. 

Since this sample is more extensive than previous studies, many effects were seen for 
the first time. Among the more remarkable was the non-detection of [C II] line in two 
normal galaxies indicating low L[cii]/Lfir, down to 2 x 10^"' (3-0" upper limits). Lower 
ratios have been seen in association with Galactic H II regions, but such low ratios were 
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unexpected in normal galaxies where the emission is from a mixture of sources: H II regions, 
WIM, and PDRs. Previous observations of galaxies had found L[cii]/Lfir to vary little, 
between 10~^ and 10~^. This result may have been biased because of the low sensitivity 
of the previous surveys and the selection criteria for the galaxies. DIRBE observations 
of the Milky Way also showed a spatially constant ratio L[cii]/Lfir=3 x 10^'^, except in 
the Galactic center. In the current sample of normal galaxies we see a smooth decline in 
L[cii]/Lfir with increasing dust temperature and star- forming activity in galaxies. In a 
sample of 60 normal galaxies, this trend spans a factor of more than 50 in L[cii]/Lfir with 
[C II] deficient galaxies at the hottest and most active end. The anticorrelation between 
L[cn]/LFiR and Fi,(60 //m)/Fi^(100 /im) is the strongest, followed by L[cii]/Lfir vs. Lfir/Lb. 
The anticorrelation between L[cn]/LFiR and the IR-luminosity is the weakest and may be 
a secondary correlation. This is good news for searches for distant startbursts. The most 
luminous galaxies are not necessarily [C II] deficient: warm galaxies are more likely to show 
a lower Lpii]/LFiR, whereas cool luminous galaxies can be detected. 

There have been many explanations for the variations in L[cii]/Lfir , but on examining 
the evidence we favor the scenario where L[cii]/Lfir decreases as the heating of the gas 
becomes less efficient in the high Gq/u regime due to charging of dust grains. As the grains 
become positively charged the efficiency of photoelectric ejection decreases. The various 
lines of evidence supporting this hypothesis are: 

(1) L[cii]/Lfir deficiency is seen in the more actively star-forming galaxies with warmer 
FIR colors. 

(2) Examination of heating and cooling balance in Arp 220 does not show any other 
lines that could be cooling the neutral ISM instead of [C II]. In our own sample, the 
inclusion of [O I] does not change the low ratio of gas to grain heating as indicated by 
(L[cii] + L[oi])/Lfir. So far we have failed to identify other channels, which leads us to 
believe that heating efficiency is low in galaxies where where L[cii]/Lfir is low. 

(3) FDR models which include grain charging and the photoelectric effect successfully 
reproduce the trends of decreasing (Ljcn] + L[oi])/Lfir and L[cii]/Lfir with 
F^(60 //m)/F^(100 //m) (Figure 8). 

The [C II] line flux shows better correlation (and less scatter) with mid-IR flux from 
galaxies which is dominated by aromatic features and emission from transiently heated 
small grains (Helou et al. 2001). This indicates that the heating of gas is dominated by 
such grains as expected by some theoretical considerations (e.g. Bakes & Tielens 1994). 

A somewhat less dramatic decrease in L[cii]/Lfir is seen in early type galaxies where 
this decrease is due to softer radiation fields. Pierini et al. (1999) also see a decrease in 
L[cii]/Lfir with decreasing star formation in quiescent galaxies. In early-type galaxies with 
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little or no star-formation, L[cii] is lower because UV photons are needed to heat the gas by 
photoelectrons, whereas both optical and UV photons can heat the dust. This scenario is 
corroborated by high ratios of Lco/L[cii] iii some early- type galaxies (Malhotra et al 2000). 
At the other morphological extreme, irregular galaxies show high L[oiii]/Lfir (Figure 6) and 
much of the [C II] seen in these galaxies may be originating in the ionized regions (Hunter et 
al. 2001). Apart from these effects, we do not see much dependence on galaxy morphology. 

Because carbon has an ionization potential lower than hydrogen, [C II] emission can 
arise in both ionized and predominantly neutral media. The observed ratio L[oi]/L[cii] is 
lower than the ratio in PDR models, suggesting a possible contribution of [C II] emission 
from ionized gas. In addition, the ratio of [C II] and [Nil] fluxes trace each other reasonably 
well in the spatially resolved data in the Milky Way (Bennett et al. 1994). These two 
facts suggest that a fair fraction of [C II] arises in ionized regions. It is also significant 
that L[oi]/L[cii] shows a remarkably tight correlation with Fi,(60 //m)/Fj,(100 /im) which is 
as expected from PDR models, suggesting that a significant fraction of [C II] arises from 
PDRs. Thus the PDRs and diffuse ionized gas are associated with each other. 

We estimate a scaling between [Nil] (122 //m) and [C II] (158 //m) from diffuse and 
dense ionized gas and use that scaling to estimate the [C II] luminosity from neutral regions 
L[cii]^. The main uncertainty in this calculation comes from the relative abundances of C 
and N. Additional uncertainty is contributed by the unconstrained fraction of the [C II] and 
[Nil] emission arising from dense vs. diffuse ionized medium, since the scaling of Ljciq/LjNn] 
is different for the two. 

By comparing the observed ratios of the [O l](63/xm) and the [CII]c lines from PDRs 
and the ratio of line to the total FIR continuum from dust with a grid of PDR models 
(Kaufman et al. 1999), we derive the average Gq to be in the range 10^ — 10'^'^, and gas 
densities n in the range 10^ — lO'^'^cm^^ for this sample of galaxies. The FUV flux Gq 
and gas densities n correlate with each other and Go increases roughly proportional to n^'"^ 
over about two orders of magnitude. We also derive Go and n from the observed values 
of L[oi]/Lfir vs Fi,(60 /im)/Fi,(100 //m). Averaged over the sample there is no systematic 
difference between Gq and n derived by the two methods. Half the galaxies show better 
than a factor of two agreement in Go and n derived from the two methods. 

The correlation between Go and n is explained by assuming that a significant portion 
of the PDR emission comes from PDRs surrounding expanding H II regions. The range in 
values of Gq and n derived for this sample suggests that the GMCs in which OB stars are 
born have different average properties from galaxy to galaxy. We suggest that the galaxies 
with higher Go and n have more massive GMCs. 
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A. Sample selection 

This sample was constructed to study star-formation and the ISM in normal galaxies. 
Here normal galaxies are deined to have the energy production in the galaxy come from 
star-formation and not active nuclei. 

The distant galaxies were selected to span a range of galaxy properties (Dale et al. 
2001): 

(1) Morphology: The morphology of the galaxies in this sample ranges from Irr through E 
(Table 5). The galaxies are uniformly distributed across the range of morphological types 
(see Dale et al. 2000). 

(2) Far-infrared luminosity: The FIR luminosity of galaxies in this sample span the range 
Log(LFiR/ Lq)=7.7-11.2. Since this is a study of normal galaxies, we avoided ultra-luminous 
galaxies which might harbor hidden AGNs. 

(3) F^(60 /im)/F^(100 /xm): The ratio of FIR fluxes in the IRAS fllters at 60 and 100 //m 
indicates the average dust temperature in the galaxies. The galaxies in the sample cover 
the range F^(60 //m)/F^(100 //m) = 0.3 - 1.37. 

(4) Lfir/Lb: The FIR to B-band flux ratio is a measure of star-formation activity and of 
dust optical depth. For low values of Lfir/Lb one can approximate Lfir as surrogate of 
extincted UV and hence current star-formation, while Lb represents relatively old stars. 
High values of Lfir/Lb indicate high star- formation rates as well as high optical depth in 
dust. 

Care was taken during sample selection to sample the parameter space in these 
four parameters: morphology, IR luminosity, F,^(60 ;um)/Fj,(100 fim) and Lfir/Lb, 
as much as possible. For a randomly selected sample, the quantities IR luminosity. 
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F,^(60 //m)/Fi,(100 //m) and Lfir/Lb, correlate with each other. We have tried to reduce 
these correlations by picking galaxies that sample the parameter space as uniformly as 
possible. Still, correlations between IR luminosity, Fj,(60 iim)/F^{100 /im) and Lfir/Lb, 
remain and are illustrated in Figure 14. 
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Fig. 13. — These two figures show how the galaxies discussed in this paper span the parameter 
space along the parameters F^(60 /im)/F,y(100 /xm), Lfir/Lb and infrared luminosity. In 
figure 13(a), a correlation between F,^(60 /im)/F;^(100 ^m.) and Lfir/Lb is seen. This 
correlation is 4.5 a significant. Figure 13(b) also shows in triangles the luminous and 
ultraluminous galaxies from the sample of Luhman et al. (1998). The correlation between 
Fi,(60 //m)/F,^(100 //m) and Lpin is 4.1(7 significant. 
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C. Physical quantities derived from compEirison with models 
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Table 6: Line fluxes for the distant sample (in units of 10 ^^W/m^) 
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< 0.029 




MRK0331 
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Table 7: Average physical parameters derived by comparing FIR line fluxes with PDR models 
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